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Abstract. Dark energy can be studied by its influence on the expansion of the 
Universe as well as on the growth history of the large-scale structure. In this paper, 
we follow the growth of the cosmic density field in early dark energy cosmologies 
by combining observations of the primary CMB temperature and polarisation power 
spectra at high redshift, of the CMB lensing deflection field at intermediate redshift 
and of weak cosmic shear at low redshifts for constraining the allowed amount of 
early dark energy. We present these forecasts using the Fisher-matrix formalism and 
consider the combination of Planck-data with the weak lensing survey of Euclid. We 
find that combining these data sets gives powerful constraints on early dark energy and 
is able to break degeneracies in the parameter set inherent to the various observational 
channels. The derived statistical ler-bound on the early dark energy density parameter 
is cr(fij) = 0.0022 which suggests that early dark energy models can be well examined in 
our approach. In addition, we derive the dark energy figure of merit for the considered 
dark energy parameterisation and comment on the applicability of the growth index 
to early dark energy cosmologies. 



PACS numbers: 98.80.-k, 95.36,+x 
1. Introduction 

Cosmological models with a dark energy fluid with an equation of state parameter w 
close to —1 are favoured by combining recent Cosmic Microwave Background (CMB), 
supernova and Baryon Acoustic Oscillations (BAO) data, suggesting that the Hubble- 
expansion accelerates in the current cosmic epoch. However, the behaviour of this dark 
energy fluid at early times is still an open question; if its equation of state remains close 
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to —1 at early times, such as for a cosmological constant, its density evolves slowly and 
was neglible relative to the matter and radiation fluids. In this case, the dark energy is 
unlikely to play a role in the physics of the early universe. 

However, in dynamic models of dark energy, such as a scalar field, this may be 
different and dark energy could have a non-negligible influence on earlier stages of the 
growth history. It has been shown that scalar field models exist with global attractor 
solutions which sub-dominantly "track" the dominant component of the cosmological 
fluid [H El El H]. In such models, the fraction of dark energy fluid to the critical 
density, Q d , was more or less constant at early times. If Q d was sufficiently large, 
its presence could have had an observable impact on the probes of the early universe, 
such as nucleosynthesis, CMB physics and the growth rate of structure. These models 
are collectively known as early dark energy (EDE) models. 

EDE models may help to resolve some apparent discrepencies in the amplitude of 
density fluctuations, usually parameterised by a 8 , the r.m.s. of the density field on scales 
of 8 h~ 1 Mpc. If the excess in the CMB power spectrum at high multipoles measured 
with the Cosmic Background Imager (CBI) [5] is explained by the Sunyaev-Zel'dovich 
(SZ) effect of unresolved clusters of galaxies one needs a value of a$ close to one, which is 
significantly higher than currently estimated from combining CMB and lensing data [6]. 
Secondly, the observed strong lensing cross section is predicted by ACDM or traditional 
dark energy models to be too small by an order of magnitude [3, El El [TO] and one 
would be forced to increase a$, again to values close to unity [II] . Analytical studies 
of spherical collapse in early dark energy models [121 03] indicate that these models 
exhibit decreased threshold redshifts for the formation of objects, which would enhance 
the SZ-signal, although recent studies question this result [T4l [i~5] . 

By combining the currently measured CMB temperature anisotropy power 
spectrum, supernova data and large-scale clustering from the Sloan Digital Sky Survey 
(SDSS), Doran et al [16] derive upper bounds on the early dark energy parameter Q d of 
0.04 (2a), in carrying out a simultaneous 8-parameter fit while employing a weak prior 
on spatial flatness and by imposing a lower bound of —1 on the dark energy equation 
of state parameter w. Linder et al [T7] confirm the relevance of early dark energy on 
the CMB power spectrum if Q d is larger than ~ 0.03, and emphasise that ignoring 
early dark energy can sever ly bias the determination of w from BAO. To improve these 
constraints, Das & Spergel [18] propose a cross-correlation of the CMB and large-scale 
structure lensing signal, from which they derive a geometrical quantity consisting of the 
ratio of angular diameter distances. Although this quantity is well suited to investigate 
dark energy models, it is rather insensitive to the value of Q d . Xia & Viel [19] use the 
latest CMB, BAO and SNIa with a prior on H Q from the Hubble Space Telescope to 
constrain the amount of early dark energy at the last scattering surface to be less than 
2% and in the structure formation era, around matter radiation equality, to be of the 
order of 6%. Adding Gamma Ray Burst and Lya forest data can tighten the constraints 
at the last scattering surface by about an order of magnitude. 

In this paper, we consider the combination of an observation of the primary CMB 
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fluctuations, the CMB lensing signal and weak cosmic shear for following the growth 
history in early dark energy models at high, intermediate and low redshift. Cosmic 
structure growth on the relevant scales is close to linear, where the theory is able to 
provide reliable predictions. The motivation of this investigation was to follow the 
cosmic growth history by combining high-precision probes at recombination redshifts, 
at intermediate redshifts of z ~ 4 by CMB lensing and at low redshifts around unity 
from weak cosmic shear. As experimental data for the microwave sky we consider the 
Planck survey, and compare to the corresponding constraints from the CMBPol satellite 
concept [20l |2T] . For the weak lensing power spectra we use the characteristics for a 
satellite mission like Euclid |22j . 

While this work was in its final stages, the work of de Putter et al [23J appeared 
which similarly uses the lensed CMB temperature and polarisation power spectra to 
constrain early dark energy at high redshifts, but combine this data with current 
supernova samples rather than future weak lensing data. They also investigate the 
degeneracies between the neutrino mass and the amount of early dark energy, because 
of their opposite influence on the dark matter power spectrum. Our predictions for the 
EDE constraints coming from lensed CMB data alone appear comparable, given the 
differences in our assumptions about the experiments. 

This article is organised as follows: the parameterisation of early dark energy models 
used in this paper is introduced in section [2] followed by a discussion of the large-scale 
structure growth. We introduce the modelling of the three observational channels in 
section [3j Constraints are presented in section H] and a summary in section [5] concludes 
the article. 



2. Early dark energy and structure growth 

2.1. Parameterisation of dark energy 

For concreteness, we use a simple phenomenological parameterisation of for a dark 
energy fluid with an early contribution proposed by Doran & Robbers 



Here is the fractional energy density today while represents early contribution. 
The EDE energy density will sub-dominantly track the evolution of the dominant 
component. The corresponding equation of state evolves from w(a -C 1) = 1/3 during 
radiation era to equality and from w(a > a e q) = during matter domination to 
w(a = 1) = w today, as illustrated by figure [TJ In true scalar field EDE models, the 
fraction of dark energy is typically proportional to (1 + wa), where Wd is the equation 
of state of the dominant component; thus, it is reduced by 25% as the universe goes 
from the radiation dominated regime to the matter dominated regime. Also, the precise 
transition of the dark energy to its late time behaviour will depend on the details of 
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Figure 1. The background evolution of the relativistic species (red), baryonic and 
dark matter (blue), and dark energy (black) are compared in the WMAP+BAO+SN 
best fit wCDM model (solid) and the EDE parameterisation (dashed) with = 0.04. 
In the lower panel we plot the evolution of the equation of state parameter w(a) 
corresponding to the two considered models. 



the scalar potential. However, while only approximate, this parameterisation should 
capture most of the essential EDE physics. 

Throughout this work we assume a flat cold dark matter model with adiabatic and 
Gaussian initial conditions. To isolate the EDE effects, we compare these models to 
those with a dark energy fluid with constant equation of state (wCDM). We assume the 
following parameters with the WMAP+BAO+SN best fit values from Komatsu et al 

m 

p = | ns = 0.96, r = 0.084, co m = tt° m h 2 = 0.1358, uj h = tt° h h 2 = 0.02267, 

tt m = 0.2732, cr 8 = 0.8049, w = -0.992, fi| = 0.03} , (2) 

while the wCDM model is identical but with Q 1 ^ = 0. Here, n s is the spectral index of 
the primordial power spectrum, r the optical depth of reionisation, u m and are the 
physical matter and baryon energy fractions, respectively, f2 m is the current fractional 
energy density of matter, a 8 is the normalisation of the matter power spectrum, w is 
the current value of the equation of state of dark energy, and is the amount of EDE 
as defined in equation ([1]). For simplicity, the helium fraction is kept fixed at Yue — 0.24. 
From these parameters we derive the Hubble constant H = 70.5 km/s/Mpc, the current 
dark energy fraction = 0.7268, the amplitude of the primordial power spectrum 
A s = 3.22 x 10~ 9 , the reionisation redshift z rcion = 10.3, and the age of the Universe 
t = 13.6 Gyr. 
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Because the equation of state of dark energy is time dependent we always consider 
perturbations in the dark energy. For an adiabatic sound speed of c 2 = 1 the fluid model 
with the above EDE parameterisation yields equivalent results as tracking scalar field 
models. In case of clustering dark energy we find that c 2 s < 1 introduces a tilt in the 
matter power spectrum in presence of EDE. The effect of changing c 2 strongly depends 
on fi|; for simplicity, here we only consider the c 2 — 1 case. 



2.2. The growth of fluctuations 

An important indicator of the presence of dark energy at any epoch is through its 
effect on the growth of fluctuations in dark matter. At late times, it is well known 
that the acceleration of the scale factor due to dark energy makes it more difficult for 
structures to collapse. However, even at early times, when the dark energy density 
scales identically to the matter density and there is no acceleration, dark energy still 
suppresses the matter growth because it does not cluster on small scales. 

On sub-horizon scales, the linear growth of structure is usually described by the 
scale- invariant dark matter perturbation normalised today, D + (a) = 5dm{a, ^)/<5dm(a- = 
1, k). It evolves according to the growth equation 



d 2 D+ 1 



da 2 a 



dlntf 

dm a 



' W+ 3 rfi m (a) J D + = 0. (3) 



da 2a 2 

The second term describes the damping of structure growth due to the expansion of 
the universe. Early dark energy increases the amount of Hubble damping without 
contributing to the density perturbation sourcing the growth, so the growth rate is 
slowed. In a standard cold dark matter model without dark energy D + oc a, but in the 
presence of early dark energy, D + oc a 1_3n d/ 5 during matter domination [26l [27]. 

This treatment is only approximate, as it does not apply on superhorizon scales 
and ignores the radiation dominated regime. In practice we use a modified version 
of the CAMB coddj] to evolve the growth numerically. In figure [2] we compare the 
evolution of the dark matter perturbation calculated using CAMB in EDE models with 
wCDM. Normalising the growth factor to its present value, we see that at redshifts 
accessible to lensing tomography, z < 3, the effect of typical EDE models on the growth 
of perturbations is at the percent level. However, at higher redshifts the presence of 
EDE is more apparent, and the integrated effect is significant. 



2. 3. The growth index parameterisation 

At late times, e.g. z < 3, the growth can be safely described by the growth equation 
([3]). The solution at these redshifts is well approximated using the well known growth 
rate formalism, or "7-index," and these approximate solutions can be a useful tool when 

I www. camb . inf o [28 
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Figure 2. Growth of matter perturbations compared for the WMAP+BAO+SN best 
fit wCDM model (blue, solid) and for the same model with additional EDE with 
= 0.04 (red, dashed) and with f2| = 0.02 (magenta, dash-dotted). We observe that 
at redshifts z < 5 (a > 0.25) the effect of EDE on the growth of perturbations is at 
the percent level. 



considering late universe observations such as cosmic shear. Specifically, in the standard 
ACDM scenario, D + (a)/a can be approximated very well through: 

OnK) 7 - 1 



D + (a)/a = exp 



-da' 



(4) 



where 7 = 0.545. When the dark energy is not a cosmological constant, this picture is 
modified somewhat. A general dark energy equation of state parameter w(a) will lead 
to a different expansion rate, and so to a different Hubble drag, leading to a shift in 7 



More generally, 7 may be considered to be a function of the scale factor or other 
cosmological parameters. Wang & Steinhardt |3T] find 

3[l-w(a)] 



7(a) 

' v ; 5-6ii;(a) 
to provide a good fit where w(a) is the effective equation of state given by: 



w(a\ 



Here, f2 m (a) is 



fim(a) 



3 f2 m (a) [1 — tt m (a 



1 - m (1 



-3uiq x 



(5) 



(6) 



(7) 
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1.4 





Figure 3. In the top panel we plot the growth factor as a function of scale factor a 
for wCDM and EDE with fl e d = 0.03, comparing the numerical solution of the growth 
equation with the approximation using the growth index 7 given in equations ([3]) and 
([5]). In the bottom part we plot the relative differences and conclude that using the 
growth index is accurate at the 10~ 3 level. 



and the prime denotes the derivative with respect to the scale factor. (While this 
expression for the matter density ignores the presence of radiation, it should be accurate 
for the low redshifts where this approximate solution is used.) For a constant equation 
of state (i.e. = 0), the growth index defined in this way is independent of the matter 
density. For more general dark energy models, and specifically when EDE is present, 7 
weakly depends on the present matter density and the scale factor. 

Linder [32] has recently suggested an alternative approach for parameterising the 
approximate solution in EDE models, introducing an overall late-time calibration factor, 
to account for the modified growth at early times due to EDE: 

7 ( w ) = 0.55 + 0.05 ^1 + w + ^nij , (8) 
<fc =1 -4.4^. (9) 

We find this approach slightly more accurate in fitting the numerical solution of the 
growth equation. In figure [3] we show the relative difference of the growth rate of dark 
matter perturbations obtained from the numerical solution of the growth equation and 
the approximation using equation (j3J) with 7 given (jSJ). The accuracy is at the 10~ 3 
level and therefore good enough for our purposes. 

Also shown in figure H] is the growth as evaluated by the CAMB code, which takes 
into account the presence of radiation. At late times the agreement is good, but the 
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scale factor a 

Figure 4. Structure growth calculated by CAMB compared with the growth equation 
in a pure matter-dark energy universe (without radiation). The error introduced by 
not including radiation is at the percent level for a < 0.1 — 0.15 (or z > 6 — 9) depending 
on the cosmological model. 

early growth is significantly altered. Beyond a redshift of 6 — 9, the radiation corrections 
become important at the percent level. 

2. 4- Impact on the density power spectrum 

The effect of EDE can be understood as a suppression of the growth of structure stronger 
than in ACDM which depends on the epoch when the scales enter the horizon. For 
modes that enter the horizon before matter-radiation equality, k > k eq , all the modes 
experience the same suppression in growth. For larger modes entering later, k < fc eq , 
the suppression does not start until the mode enters the horizon which leads to a tilt 
in the power spectrum at large scales, see Caldwell et al [33] . In figure [5] we plot the 
matter power spectrum for the fiducial model and compare it to the EDE model with 
the same cosmological parameters plus VL e A = 0.04. 

Although the shape of the matter power spectrum can be measured by galaxy 
surveys, the suppression due to EDE at the accessible scales is degenerate with the 
galaxy bias and the normalisation erg, as can be seen in figure El thus, EDE cannot be 
constrained with these observations alone. At the same time, nonzero f2| introduces a 
tilt in the power spectrum at large scales, k < k eq , and is therefore degenerate with the 
primordial spectral index n s . Thus the combination of information on the growth of 
structure at different redshifts and at different scales is needed to determine a possible 
EDE component. 
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3. CMB, CMB lensing and cosmic shear 

3.1. Information from the CMB 

The CMB is most sensitive to the physics of the early Universe, and observations of 
the CMB temperature and polarisation anisotropies are primarily useful for fixing the 
parameters, n s , r, u m , c<j b , as well as the primordial amplitude A s . In the usual wCDM 
case, the sensitivity to the dark energy parameters (wq and Q m , which gives the present 
dark energy density) comes entirely through the observed angular scale of the sound 
horizon, since late dark energy changes the expansion and thus the distance to the last 
scattering surface. 

In EDE models, there is significant dark energy present even at last scattering. The 
primary effect of this is to change the time of last scattering, and thus the actual sound 
horizon scale (in addition to affecting the distance to the last scattering surface) [35| l36]. 
Since only one observable is affected, the angular scale of the sound horizon, the various 
dark energy parameters wo, fl m } are largely degenerate using CMB data alone. 
(Note that the apparent sensitivity to f2| seen in the figures below is a result of our 
normalisation to erg, which is not probed directly by the CMB.) 

Information from CMB lensing and/or from other observational channels, like weak 
cosmic shear, are needed to break the degeneracies and measure these late Universe 
parameters accurately. These observations provide a direct probe of the the present 
amplitude of perturbations, related to <r 8 , and by combining these with the CMB 




Figure 5. The matter power spectrum is plotted for the WMAP+BAO+SN best 
fit wCDM model (blue, solid) and for the same model with EDE (red, dashed). We 
also plot the non-linear correction using the Halofit fitting formula based on N-body 
simulations by Smith et al [31] for the best fit model (black, thin solid) and the EDE 
model (magenta, thin dashed). In the left panel as was fixed to the fiducial value 
by adjusting A s accordingly when EDE was added. In the right panel A s was held 
constant while a$ = 0.6718 for tt e d = 0.04. 
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constraints on A s one can constrain the amount of growth between last scattering and 
the present. This provides a large lever arm to probe the suppression of growth caused 
by the presence of EDE. Lensing measurements have the advantage that they are not 
subject to the bias issues which plague galaxy clustering observations. 

3.2. Information from CMB lensing 

CMB lensing is the gravitational deflection of CMB photons by the large scale structure 
between the last scattering surface and the observer. This secondary signal in the CMB 
has several statistical effects on the CMB anisotropies which are detectable. In this 
context we are interested in the changes to the power spectra and the generated 5-mode 
polarisation, which can be used to extract cosmo logical information. See e.g. Lewis & 
Challinor [37] for a comprehensive review of CMB lensing and Smith et al [38] and 
Perotto et al [39] specifically on the estimation of cosmological parameters from the 
lensed CMB spectra. 

The total deflection angle of a photon on its way to the observer is of the order of 
2 arcminutes and therefore we might expect lensing to become an important effect at 
I > 3000. Unfortunately, on these angular scales, many foregrounds exist which could be 
problematic for the interpretation of the lensing signal. However, the deflection angles 
are correlated over the sky by an angle of ~ 2 degrees, which is about the scale of the 
primary acoustic peaks in the CMB. This means that degree sized cold and hot spots 
appear after lensing larger or smaller by ~ 2', or ~ 3%. On average, lensing changes 
the statistics of the size distribution qualitatively, leading to a ~ 3% broadening of the 
acoustic peaks and dominates the spectrum on arcminute scales. 

The lensing effect is described by re-mapping the CMB anisotropies by the deflection 
angle ex. such that the lensed temperature in direction n, @(n), is given by the unlensed 
temperature in the deflected direction, 6(n) = 0(n + a). The tilde denotes lensed 
quantities. At linear order perturbation theory the deflection angle can be calculated 
as the gradient of a projected potential, ex = V^. The so-called lensing potential ip is a 
weighted integral of the gravitational potential along the line of sight and depends via 
the lensing efficiency function on cosmology. 

The interesting property of the lensing potential is that it has contributions from 
the structure out to quite high redshift. And in fact it is only weakly sensitive to 
late time and non-linear evolution. This characteristic makes CMB lensing a powerful 
tool for investigating EDE. In figure [6] we compare the primary CMB temperature and 
polarisation power spectra and the power spectrum of the lensing potential for the 
fiducial model with those in case we add a small amount of EDE. Because we normalise 
to o" 8 today, overall, EDE enhances the amplitudes of the CMB power spectra as the 
perturbations suffer a larger Hubble drag during their evolution. The lensing potential 
power spectrum is a weighted integral of the matter power spectrum and therefore also 
experiences a scale dependent enhancement at large scales (compare left panel of figure 
[5] and lower right panel of figure [6]) . 
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Figure 6. The CMB power spectra for the WMAP+BAO+SN best fit wCDM model 
(blue, solid) are compared with those in case of EDE with = 0.04 (red, dashed). 
Because we normalise to <jg today, the most apparent effect of EDE is an enhancement 
of power at the time of recombination. 

The lensing potential, however, is not directly observable and thus, in an 
experiment, has to be estimated via reconstruction techniques using other observables. 
Lewis & Challinor [37J give an overview over the different possible techniques for the 
extraction of the lensing signal. All lensed auto- and cross-correlation power spectra 
of the temperature and the polarisation can be used to construct quadratic estimators 
for the lensing potential power spectrum as proposed by Zaldarriaga & Seljak [40J and 
Guzik et al [41J. These are then combined to a minimum variance quadratic estimator 
as introduced by Hu & Okamoto [12] using a flat sky approximation. Okamoto & 
Hu [43J construct it on the full sky which we implement to estimate the noise on the 
reconstructed lensing power spectrum. 

3.2.1. CMB Fisher matrix calculations To explore the power of CMB lensing to 
measure the growth history at intermediate redshifts and therefore to break the 
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degeneracies between the late Universe parameters, we perform a Fisher matrix forecast 
for the upcoming Planck mission and later compare it with the futuristic proposal 
CMBPol. The Fisher matrix formalism is a straightforward and computationally 
fast method for predicting the precision with which future observations can measure 
cosmological parameters in a given model. The main assumption is that the parameter 
likelihood can be approximated by a Gaussian close to the maximum. 

The Fisher matrix for a CMB experiment including lensing information is 
comprehensively presented by Perotto et al [39] and the results are compared to those 
resulting from sampling the likelihood by Monte Carlo Markov-chain methods. The 
errors derived from the Fisher formalism are in very good agreement with the Monte 
Carlo results for most model parameters. For strongly degenerate parameters and where 
parameter space is bounded from below or above the Fisher matrix is still found to yield 
a satisfactory estimate for an upper bound on the error. 

The likelihood function £(data|p) is approximated as a Gaussian function of the 
model parameters p. The likelihood is Taylor expanded about its peak at p = p°. Then 
the Fisher information matrix is defined as the relevant second order term 



According to the Cramer-Rao minimum variance bound, the variance of the parameter 
p a , if marginalised over all other parameters, is given by: 



For a CMB experiment the theoretical predictions are given in terms of the power 
spectra Cf p ' where PP' £ {TT, EE, TE, ipvp, ipT}. Then the Fisher matrix can be 
expressed in terms of the derivatives of the theoretical power spectra with respect to 
the model parameters as 

1=2 PP',QQ' ^ a ^ p 

Covi is the covariance matrix of the power spectra at multipole £. It is constructed from 
the observed power spectra, including all sources of variance to the fields, see Perotto 
et al [39] for details. For the CMB temperature and polarisation we add homogeneous 
instrumental noise and beam as given by Knox [H] 



where A P is the detector noise for P £ {T, E, B} and ^fwhm is the full width at half 
maximum of the beam. In the table [T] we give the characteristics of the different CMB 
experiments considered. Generally we take the maximal multipole to be £ max = 2250 
for both CMB experiments. Planck is noise dominated in all bands already for smaller 
multipoles and our errors from CMBPol will be an upper bound. 

In figure [7J we plot the different theoretical power spectra for the fiducial model 
and compare them with the experimental noise. In case of the lensing potential power 




(10) 




(11) 




(13) 
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Figure 7. The CMB power spectra for the WMAP+BAO+SN best fit wCDM model 
are plotted against the expected noise levels for Planck and CMBPol. The lower panels 
show the cumulated signal-to-noise as a function of the highest multipole. 



spectrum the noise is estimated as described above using the noisy temperature and 
polarisation power spectra given here. We also plot the cumulated signal-to-noise which 
we define as 

s/m s E (2f 2 +1) c ^V <") 

£'=2 

To calculate the noise power spectra and the Fisher matrix for the Planck and 
CMBPol, we extended the Boltzmann code CAMB to include the EDE parameterisation 
defined in equation ([I]) and implemented the lensing reconstruction on the full sky from 
Okamoto & Hu p]. 

In figure [8] we plot, first, the derivatives of the primary temperature and E-mode 
power spectra with respect to the dark energy parameters wo and f^. Second, we plot 
the derivative of the power spectrum of the lensing potential with respect to all model 
parameters. 

We checked the numerical stability of the derivatives (finite differencing) by varying 





Planck 


CMBPol 


temperature noise Ay 
polarisation noise Ae,b 
beam $fwhm 
fraction of sky / sky 


28 /zK arcmin 
57 /iKarcmin 
7 arcmin 
0.65 


1 /iK arcmin 
1.4 /iKarcmin 
3 arcmin 
0.65 



Table 1. Assumptions on the characteristics of the different CMB experiments 
considered. Generally we take the maximal multipole to be ^ max = 2250 for both 
CMB experiments. 
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Figure 8. Left panel: logarithmic derivatives of the temperature, -E-mode, and lensing 
potential power spectra with respect to (dashed) compared to those with respect to 
wq (solid). Right panel: logarithmic derivatives of the lensing potential power spectrum 
with respect to all model parameters. 



the step size, comparing the usual double sided derivative with the one sided, and 
by using a 5-point derivative. No significant change in the results was observed. 
Furthermore, we calculated the results at different fiducial values for f2| and found 
this to be up to a 10% effect on the marginalised errors. 

3.3. Information from cosmic shear 

We derive the sensitivity of typical next-generation tomographic weak lensing surveys to 
cosmological parameters following the recent literature [45J, HfH H3 EHJ [49] . In particular, 
we study a Euclid type survey as laid out by Amara & Refregier [50], which dealt with 
wCDM, and for modified gravity theories we refer to Amendola et al [5"T] . 

Let us briefly recall the main equations for weak lensing studies. The weak lensing 
convergence power spectrum (which in the linear regime is identical to the ellipticity 
power spectrum) is a linear function of the matter power spectrum convolved with the 
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lensing properties of space; and it can be written as 



/OO IT / tt n \ 

J^WMWMP* [Pun (^J, 



dz (15) 



where P n i[Pu n (k, z)\ is the non-linear matter power spectrum at redshift z obtained by 
correcting the linear power spectrum Pn n (k,z). Here, the subscripts i and j, refer to 
bins in redshift space of the lensed galaxies, which are most commonly separated using 
redshifts obtained photometrically (photo-z's.) In flat space we have: 

W t (z) = ^a a F,(z)(l + z) , (16) 
JZi r{0, z s ) 

POO 

n t {z) = Di(z)f / D^dz', (18) 
Jo 

^-Cm" (19) 

where Di(z) is the radial distribution function of galaxies in the i-th redshift bin. We 
assume an overall radial distribution 

D(z) oc z 2 exp [-(z/z ) 3/2 ] . (20) 

The distributions Di are obtained by binning the overall distribution and convolving 
with the photo-z distribution function. 

3.3.1. Cosmic shear Fisher matrix calculations The Fisher matrix for weak lensing is 
given by 

sr (2£ + l)A£ dP ij t dP km x 
F aP - / sky 2. j hk ^v~i { U [ } 

and the covariance matrix is defined as 

C jk = P 3k + 5 3k ( 1 l t )nj 1 (22) 

where 7 int is the r.m.s. intrinsic shear (we assume (7^) = 0.22 [50J) and 

nj = 3600 d (^pj hj (23) 

is the number of galaxies per steradian belonging to the z-th bin, d being the number of 
galaxies per square arcminute and hi the fraction of sources belonging to the z-th bin. 

In the following we consider a satellite mission like Euclid with the characteristics 
given by Refregier [22]: sky fraction / s k y = 1/2, mean redshift of the galaxy distribution 
-Zmean = 0.9, and number of sources per arcmin 2 of d — 40. We assume that the photo- z 
error obeys a normal distribution with variance a z = 0.05(1 + z). We choose to bin the 
distribution out to z = 3 into five equal-galaxy-number bins. 

In figure [9] we plot the convergence power spectrum for the fifth redshift bin, P^(£) 
for EDE fiducial model, with the noise due to the intrinsic ellipticity and for comparison 
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Convergence Power Spectrum (5-5 bins) 




Figure 9. The convergence power spectrum for the fifth redshift bin, Pss^), is shown 
for the WMAP+BAO+SN best fit wCDM model (blue, dashed) compared with the 
EDE fiducial model (red, solid). The lower panel shows the relative difference between 
the two. The error bars represent the noise errors on the EDE convergence power 
spectrum in logarithmically spaced bins. 

we also plot the wCDM spectrum. We see that the EDE spectrum and the wCDM 
spectrum are almost undistinguishable. The derivatives of the same convergence power 
spectrum with respect to the model parameters are shown in figure [H]] 

For the linear matter power spectrum we adopt the fit by Eisenstein Sz Hu [53] 
(with no massive neutrinos and also neglecting any change of the shape of the spectrum 
for small deviations around w = —1). For the redshifts of interest, we verified that the 
Eisenstein & Hu fit is accurate enough at the scales to which cosmic shear is sensitive. 
For the non-linear correction we use the fitting function to the halo model from Smith 
et al [31]. We consider the range 10 < t < 10000 since we find that both smaller and 
larger multipoles do not contribute significantly. 

For all results presented in the next section we assumed that the systematic 
errors can be neglected. A partial list of these errors include, on the observational 
side, shape measurement systematics and photometric redshift systematics and, on the 
theoretical side, the uncertainties concerning the accurate prediction of the small scale 
matter spectrum (due to non-linear effects and to the uncertain baryon contribution), 
the uncertainty on the source distribution, the degeneracies with other effects like 
massive neutrinos. Some of these effects have been shown to be under control, e.g. the 
photometric redshift systematics and, to some extent, the shape measurements [50J. 




4. Results 

4.1. Results for CMB and CMB lensing 

In figure [TT] we present the forecasted la error ellipses of all combinations of the 
model parameters for three cases: first, Planck without lensing information, second, 
Planck including the reconstructed lensing power spectrum, and last, CMBPol also 
including lensing information. The standard deviation on each single parameter (fully 
marginalising over all other parameters) is given in table [2] for the three cases. 

The lensing information in Planck, as expected, does not improve the constraints 
on the early Universe parameters {n s , r, uj m , c^b} significantly. On the other hand, 
it starts to break the degeneracy between the late Universe parameters {fi m , cr 8 , w } 
slightly. Finally, it improves the constraints on Q 1 ^ by a factor ~ 1.4, and so Planck will 
be able to constrain EDE at sub-percent level, cr(f2^) = 0.0037. 

CMBPol, the futuristic proposal of a dedicated satellite mission to map out the 
CMB polarisation modes, will obviously be able to constrain all parameters very well, 
as also pointed out by Smith et al [21] and de Putter et al [23]. The constraints on the 
single parameters are improved by factors ~ 1.8 — 4 compared to Planck with lensing. 
Note that this result strongly depends on the achievable noise levels and the maximal 
multipole of CMBPol. 

4-2. Results for cosmic shear 

In figure [12] we show the confidence regions for the different combinations of the late 
Universe parameters {fl m , ^s, w , Q^} . We always plot the fully marginalised confidence 
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Figure 11. Error ellipses (la) for all combinations between the late Universe 
parameters while marginalised over the other parameters. Projected constraints for 
Planck without CMB lensing are dotted (yellow), for Planck including CMB lensing 
are dashed (light green), and for CMBPol are solid (dark green) ellipses. Clearly 
already the CMB lensing information from Planck will help to break the degeneracies 
between late time parameters slightly. 



ellipses at 68%, which in two dimensions corresponds to semiaxes of length 1.51 times 
the eigenvalue. Marginalising over all other parameters, we see that errors of 0.04 for 
wq and 0.024 for f2| are achievable with a Euclid type survey. In table [2] the la errors 
on all model parameters are listed. Since the matter power spectrum does virtually not 
depend on the optical depth of reionisation a cosmic shear survey does not constrain r. 

The predicted error cr(Q^) is comparable to the recent upper limit of f2| < 0.045 
by Bean et al |54j indicating that galaxy weak lensing experiments are not strongly 
sensitive to an early dark energy component. On the other hand, the strong constraint 
on Wq from cosmic shear tomography will help to constrain dark energy in combination 
with information from the CMB as we discuss in the next section. 

In principle, one could consider a more nearer-term weak lensing survey, like 
CFHTLS Widqa or However, their sky coverage is considerably smaller and 

therefore the statistical errors significantly larger. We verified that CFHTLS would 



§ The Canada-France-Hawaii Telescope Legacy Survey: www.cfht.hawaii.edu/Science/CFHLS 
|| The Dark Energy Survey: www.darkenergysurvey.org 
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Figure 12. Confidence regions at 68% for the a Euclid type lensing tomography 
survey with z mean = 0.9 and assuming 40 galaxies per arcmin 2 . 

not and DES would only marginally improve the constraints on EDE form Planck, 
respectively. 

4-3. Combined results 

We combine the Fisher matrices for the possible CMB, CMB lensing and cosmic shear 
measurements. The redshifts which are probed by CMB lensing and cosmic shear 
tomography are quite distinct from each other such that we can safely neglect the 
covariance between the measurements. This means we simply use the sum of the two 
Fisher matrices to calculate the combined constraints. 

In table [2] we summarise the la uncertainties derived for Euclid, Planck, and 
CMBPol, and for the combinations of Planck+Euclid and CMBPol+Euclid (both 
incl. CMB lensing). The combinations clearly benefit from the information on the 
early Universe from the CMB and on the late Universe from cosmic shear. From 
Planck+Euclid we infer a constraint on EDE of a(Q^) = 0.0022 when marginalised over 
all other parameters. This improves the constraint from Planck alone by a factor 1.7. 
In figure [13] we plot the la confidence contours from Euclid, Planck, and Planck+Euclid 
for the different combinations of all parameters except r which is not constrained by 
Euclid. 
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Table 2. Forecasted errors (ler) for the model parameters from the CMB: first, 
for a Euclid type space based cosmic shear survey, second, for Planck without 
lensing information, third, for Planck with lensing information reconstructed from 
the primary power spectra, fourth and fifth, for the futuristic proposal CMBPol with 
and without lensing reconstruction, then, for the combinations of Planck+Euclid and 
CMBPol+Euclid (both incl. CMB lensing). Note that Planck will already be able to 
constrain EDE to sub percent level and the combination with Euclid improves the 
constraints on £1^ and wq to a level comparable with CMBPol. 



In comparison, CMBPol alone will be able to constrain down to 0.0013. Overall, 
Planck+Euclid does quite well compared with CMBPol. This is, first, because the 
marginalised constraints on the late Universe parameters {fl m , a 8 , w } from Euclid are 
comparable with those from CMBPol. And second, because the likelihood surfaces in 
the (wq, Q 1 ^) plane for Planck and Euclid are strongly complementary as we point out in 
figure [TU The combination CMBPol+Euclid would enhance all constraints even more. 

4-4- EDE figure of merit 

To be able to compare the leverage of the considered experiments on EDE we propose a 
figure of merit (FOM) based on the strength in constraining wq and fl^ simultaneously 
while marginalising over all other parameters. We define F _1 as the 2x2 sub-matrix of 
the inverse Fisher matrix where the rows and columns associated with the marginalised 
parameters have been removed. Then we define 

FOM= VdetF. (24) 

In table [3] the FOMs for the different experiments in consideration are listed together 
with the forecasted la constraints on the dark energy parameters. Although CMBPol 
alone (incl. lensing) does much better than Planck, the FOM of the combination 
Euclid+Planck is a factor ~ 1.4 higher than the FOM of CMBPol. This is because 
of the strong constraints from Euclid on the late-universe parameters which anchor the 
CMB and CMB lensing constraints at low redshifts. 
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Figure 13. Error ellipses (la) for all model parameters for CMB lensing from 
Planck (green, dashed), for cosmic shear from Euclid (turquoise, dotted), and for 
their combination (blue, solid). 



5. Summary & conclusions 

The focus of our investigation has been on early dark energy, its impact on the cosmic 
growth history and its observability with three structure formation probes, specifically 
the primary CMB temperature and polarisation spectrum at high redshift, the CMB 
lensing deflection field at intermediate redshift and weak cosmic shear at low redshift. 

• Early dark energy models show a scale dependent suppression of the cosmic density 
field during horizon entry and cause a tilt in the matter power spectrum, which 
is naturally degenerate with n s and erg. The Newtonian growth in early dark 
energy cosmologies is considerably higher than in standard dark energy models 
with varying equation of state. The growth function can be well approximated by 
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Figure 14. Error ellipses [la) for the dark energy parameters wq and 51^ from 
Planck (green, dashed), for cosmic shear from Euclid (dotted turquoise), for the 
combination Planck+Euclid (blue, solid), and from CMBPol alone (dark green, solid). 
The likelihood surfaces in the (u>o,f2§) plane for Planck and Euclid are strongly 
complementary which leads to a great enhancement of the constraints when combining 
the two probes. 



Experiment 


FoM 


a(w ) 




Euclid 


2143 


0.040 


0.024 


Planck 


1500 


0.13 


0.0051 


Planck lensing 


3210 


0.090 


0.0037 


CMBPol 


16075 


0.026 


0.0025 


CMBPol lensing 


35629 


0.023 


0.0013 


Euclid+ Planck 


49765 


0.016 


0.0022 


Euclid+ CMBPol 


157213 


0.0080 


0.00092 



Table 3. We compare the FOM defined in equation (|24jl and the la errors on wq 
and fij for the different experiments and the combinations of Euclid with Planck 
and CMBPol (both incl. lensing). Note that Planck+Euclid has a higher FOM than 
CMBPol alone. 

the growth index. 

• The differences in the growth function at high redshift and the modifications of 
the power spectrum on large spatial scales suggest that the CMB fluctuations and 
CMB lensing are well suited for investigating early dark energy. Weak cosmic shear 
as a low redshift probe provides a precision measurement of fl m , as and w and is 
very useful for breaking degeneracies 
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• We compute constraints on the early dark energy density and the degeneracy of Vt e d 
with other cosmological parameters from the CMB temperature and polarisation 
power spectra, with corresponding noise levels and angular resolution for the Planck 
surveyor and for the proposed CMBPol mission. The Fisher analysis suggests 
that the statistical uncertainty on Vt e d is o~(Q d ) = 0.0051 from the CMB spectra 
observed by Planck, which decreases to cr(f2|) = 0.0037 including CMB lensing. 
The respective numbers for CMBPol are a(Q e d ) = 0.0025 and a(Q e d ) = 0.0013. We 
improve these measurements and break degeneracies by combining them with the 
Euclid weak lensing survey. By effectively adding priors on Q m , a 8 and w at low 
redshifts, weak lensing data further decreases the uncertainty to c(Q d ) = 0.0022 
for Planck and o~(Q d ) = 0.00092 for CMBPol. The orientation of the degeneracy 
ellipses suggests that weak lensing is very well suited for breaking degeneracies and 
hence improving the constraints on Q d and Wq. 

• We define a figure of merit defined as the inverse volume of the parameter space 
spanned by fl d and the dark energy equation of state w contained inside the la 
isolikelihood contour. The increase of accuracy by combining CMB and lensing 
data is very illustrative: both Planck and Euclid for themselves reach values of 
about 2 — 3 x 10 3 , whereas their combination attains almost 5 x 10 4 . A similar but 
lower value is reached for the CMBPol CMB polarisation power spectra, suggesting 
that a precision determination of these two cosmological parameters is possible with 
upcoming data sets. 

• In principle, other observational channels are also sensitive to While the 
distance-redshift relation from supernovae is only expected to serve as a prior on 
the late Universe parameters, the integrated Sachs- Wolfe effect is mildly sensitive 
to the a possible early evolution of dark energy. However, this will be difficult to 
measure. More interesting will be future observations of the 21 cm radiation from 
reionisation with LOFAR and eventually with SKA. We expect a similar leverage 
on early dark energy from 21 cm lensing as from CMB lensing, see e.g. Jarvis [55] 
and Metcalf & White [56]. 
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